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Abstract 

We discuss the quasiclassical Green function method for a two-dimensional electron gas in the presence of spin-orbit coupling, 
with emphasis on the meaning of the ^-integration procedure. As an application of our approach, we demonstrate how the 
spin-Hall conductivity, in the presence of spin-flip scattering, can be easily obtained from the spin-density continuity equation. 
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1. Introduction 

The quasiclassical technique is one of the most pow- 
erful methods to tackle transport problems. Its main 
virtue relies in the fact that starting from a micro- 
scopic quantum formulation of the problem at hand it 
aims at deriving a simpler kinetic equation resembling 
the semiclassical Boltzmann one. In deriving such an 
equation some of the information at microscopic level 
is suitably incorporated in a set of parameters charac- 
terizing the physical system at macroscopic level. Since 
the first application to superconductivity, this equa- 
tion is known as the Eilenberger equation (for a review 
see for instance [1]). We have recently derived[2] such 
an equation for a two-dimensional electron gas in the 
presence of spin orbit coupling with Hamiltonian 
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where b(p) is a momentum dependent internal mag- 
netic field. In the case of Rashba spin-orbit coupling 



b = ap A e z . In Ref. [2] we adopted the standard £- 
integration procedure to arrive at the Eilenberger equa- 
tion, and, though this leads to correct results, we feel 
the need for a deeper understanding, which we provide 
in the present paper. In so doing we follow an analysis 
carried out by Shelankov[3]. Finally, we use the Eilen- 
berger equation to study the response to an external 
electric field in the presence of magnetic impurities. 



2. The quasiclassical approach 

In deriving the Eilenberger equation a key observa- 
tion is that, by subtracting from the Dyson equation 
its hermitian conjugate, one eliminates the singular- 
ity for equal space-time arguments and gets a simpler 
equation for the ^-integrated Green function 
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Here G(p, x) is the Green function in Wigner space, 
i.e. the Fourier transform of G(xi,X2) with respect to 
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the relative coordinate r = xi — X2 . The "check" indi- 
cates that the Green function is a 2 by 2 matrix in the 
Keldysh space [1]. To shed some light on the meaning 
of the ^-integration, let us consider first the space de- 
pendence of the two-point retarded Green function for 
free electrons in the absence of spin-orbit coupling 



G fl (x 1 ,x 2 ) = ^ — 



Xi - x 2 . 



(3) 



At large distances, the integral is dominated by the ex- 
trema of the exponential under the condition of con- 
stant energy. This forces the velocity to be parallel or 
antiparallel to the line connecting the two space argu- 
ments, <9 P e(p) oc r. It is then useful to consider the mo- 
mentum components parallel (py) and perpendicular 
(p±) to r. Given the presence of the pole, one can ex- 
pand the energy in powers of the two momentum com- 
ponents £(p\\,p±) = vf(p\\ — Pf) +p 2 ± /2m. In the case 
of the retarded Green function, the important region 
is that with velocity parallel to r. We then get 
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One sees how the Green function is factorized in a 
rapidly varying term ~ e lpFr / \/ppr, and a slow one, 
e i(ui/v F )r _ -pjjjg SU gg es ts to write quite generally 



G ii (x 1 ,x 2 ) = - A /— AT e ipFr ff fl (xi,x 2 ) 
V Pft 
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where g fl (xi,x 2 ) is slowly varying and G R indicates 
the free Green function. Explicitly, in the present equi- 
librium case 
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For the advanced Green function one can go through 
the same steps with the difference that the integral is 
dominated by the extremum corresponding to a veloc- 
ity antiparallel to r, so that one has the ingoing wave 
replacing the outgoing one. In the non-equilibrium case 
Shelankov has shown that 

g R { W2 ) = i-f dCe^G ii (p,x), p = pr (7) 



and furthermore that the quasiclassical Green function 
corresponds to the symmetrized expression 



ff *(fcx) = Mm - / d£cos (£) G fl (p,x) 
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when sending to zero the relative coordinate r. 

When the spin-orbit coupling is present the Green 
function becomes a matrix in spin space and the Fermi 

2 

surface splits into two branches e±(p) = §^ ± |b|. We 
always assume this splitting to be small compared to 
the Fermi energy, i.e. |b|/eF <Si 1. In the case of the 
Rashba interaction we write 



G«(x 1 ,x 2 ) = -X: v /giiV ± e- ± ^{p ± ,^(xi ; x 2 )} 

(9) 

where P± = |±)(±| is the projector relative to the 
± energy branch and the curly brackets denote the 
anticommutator. This ansatz allows us to proceed in 
Wigner space as before, while retaining the information 
on the coupling and coherence of the two bands. Eq.(9) 
is the equivalent in real space of the ansatz for the 
Green function G(p, x) used in Ref.[2]. With such an 
ansatz, Eq.(9), we obtain from Eq.(7) 

S K (x 1; x 2 ) = J2 ^ {iV ± P ± ,ff fl ( Xl ,x 2 )} . (10) 

What we have explicitly shown for the retarded com- 
ponent of the Green function can be extended to the 
advanced and Keldysh components too. Notice that 
g R and g R coincide in the absence of spin-orbit cou- 
pling, since in that case iV± = iVo- The derivation of 
the Eilenberger equation can now be done following the 
steps detailed in Ref.[2]. We do not repeat them here 
and give just the final result 

y~] {dtfjv + \\ ~^ + d P^° v ' "O'dx^j + i[b„ ■ cr,g„]) 

= -i[%§], (11) 

where g v = (1/2){P„, g}, g = g + + <?_ and both the 
momentum p„ and the internal field b„ are evaluated 
at the !/-branch of the Fermi surface. Finally, E is the 
self-energy. It is often convenient to expand g in terms 
of Pauli matrices, g = go + g • cr , to explicitly separate 
charge and spin components. Physical quantities like 
charge and spin densities and currents are related to the 
Keldysh component of g. For example the spin current 
for si, I = x, y, z is 
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U*,t) = -\*N j gjf'Cejx.t), (12) 



where 



and (...) is the angle average over the directions of p. 



3. Magnetic impurities and spin currents 

Focusing on the Rashba interaction, we study the 
effects of magnetic impurities on spin currents. In [4] 
and [5] the problem has been recently tackled via dia- 
grammatic techniques. We show how analogous results 
can be obtained in a simple and rather elegant way 
relying on eq.(ll). As it is well known, spin currents 
arising from the spin Hall effect in such a system are 
completely suppressed by the presence of non-magnetic 
scatterers. By taking the angular average of eq.(ll), 
one obtains a set of continuity equations for the vari- 
ous spin components which let one easily understand 
the origin of this cancellation. Explicitly, by assum- 
ing s-wave and non-magnetic impurities randomly dis- 
tributed in the system 



Vi(x) =J2 U6(x--Ri), 



(14) 



the self-energy in the Born approximation turns out to 
be E i = —i{g}/2r, 1/r being the momentum scatter- 
ing rate. The continuity equations for the I — x,y,z 
spin components then read 

9t<5()+a x - Ji = 2(b Ag)i. (15) 

A rather important peculiarity of the Rashba Hamil- 
tonian is that it lets one write the vector product ap- 
pearing above in terms of the various spin currents, so 
that, by choosing for example I = y, we are left with 



9t{g y ) + <9 X ■ Js = -2maJs,y. 



(16) 



Under stationary and homogeneous conditions this im- 
plies the vanishing of the J* y spin current. As soon 
as magnetic impurities are introduced in the system, 
their presence changes the self-energy and leads to the 
appearance of additional terms in Eq.(16). We assume 
the magnetic scatterers to be also isotropic and ran- 
domly distributed 

V 2 (x) = J2 B-<r<5(x-Ri), 



and, proceeding again in the Born approximation, we 
obtain the self-energy 
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£ = £ 1+ £ 2 = --L(g)-^J2 °i<SlW (18) 

Here l/r s f is the spin- flip rate. With this, and by con- 
sidering again stationary and homogeneous conditions, 
Eq.(16) becomes 



2maJ*y + -^—{9y) = 0, 



(19) 
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which in terms of the real spin current and polarization 
means 



V 



(20) 



By assuming a low concentration of magnetic impuri- 
ties, we can use in eq.(20) the value of the {/-spin polar- 
ization valid in their absence, s y — —\e\EarNo[6], E 
being the external, homogeneous electric field. We then 
get the spin Hall conductivity to first order in t/t s / 

OsH - ~ , (21) 

a results that differs from those on Refs. [4,5] . This is not 
surprising for Ref . [5] , which neglects normal impurity 
scattering and then considers the opposite limit. The 
reason why our result does not agree with the low mag- 
netic impurity-concentration limit of eq.(20) of Ref.[4] 
is not clear to us and deserves further investigation. 
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